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In this paper we consider the density, at a point z — x + iy, of critical percolation clusters that 
touch the left (Pl(z)), right (Pr(z)), or both (Plr(z)) sides of a rectangular system, with open 
boundary conditions on the top and bottom sides. While each of these quantities is nonuniversal 
and indeed vanishes in the continuum limit, the ratio C(z) = Plr(z)/ ^/ Pl(z)Pr(z)IIi x , where II^ 
is the probability of left-right crossing given by Cardy, is a universal function of z. With wired 
(fixed) boundary conditions on the left- and right-hand sides, high-precision numerical simulations 
and theoretical arguments show that C(z) goes to a constant C = 2 7/2 3" 3/4 tt 5/2 r(l/3)~ 9/2 = 

S 1.0299268 . . . for points far from the ends, and varies by no more than a few percent for all z values. 
Thus Plr(z) factorizes over the entire rectangle to very good approximation. In addition, the 
numerical observation that C(z) depends upon x but not upon y leads to an explicit expression 
for C(z) via conformal field theory for a long rectangle (semi- infinite strip). We also derive explict 
expressions for Pl(z), Pr(z), and Plr(z) in this geometry, first by assuming y-independence and 
then by a full analysis that obtains these quantities exactly with no assumption on the y behavior. 
In this geometry we obtain, in addition, the corresponding quantities in the case of open boundary 

S conditions, which allows us to calculate C(z) in the open system. We give some theoretical results 
for an arbitrary rectangle as well. Our results also enable calculation of the finite-size corrections to 
the factorization near an isolated anchor point, for the case of clusters anchored at points. Finally, 
we present numerical results for a rectangle with periodic b.c. in the horizontal direction, and find 
Q C(z) approaches a constant value Ci ~ 1.022. 
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qq I. INTRODUCTION 

o ' 

Percolation at the critical point has many well-known universal properties, including universal critical exponents, 
scaling functions, and amplitude ratios. Universality means that the properties are the same for all realizations of the 
system (with a given dimensionality) in the continuum or field-theory limit. Crossing probabilities are also universal 
and have received a great deal of attention since the work of Cardy [l[ and Langlands et al. and renewed interest 
more recently with the development of Schramm- Loewner Evolution (SLE) 01 and a new set of results for percolation 

1- 

A more detailed picture of the critical system can be obtained by examining the correlations within clusters of 
connected sites. In previous work 0, we demonstrated, by use of conformal field theory and high- precision 
simulation, certain exact and universal factorizations of higher-order correlation functions in terms of lower-order 
correlation functions for percolation clusters in two dimensions at the percolation point. In that work, the correlation 
functions involved the density of critical percolation clusters constrained to touch one or two isolated boundary points, 
or single boundary intervals. Here we extend those results by considering densities constrained to touch one or two 
distinct boundary intervals, which is a more difficult problem. 

Specifically, we consider the quantities Pl(z), Pr(z), and Plr{z), which give the density of percolation clusters at 
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a point z — x + iy that touch the left, right, or both sides of a rectangle, respectively, as well as IT/,, the probability of 
a horizontal crossing (i.e., one or more clusters that touch both left and right sides) which is given by Cardy's formula 
Pl(z), Pr(z), and Plr(z) also determine the probabilities that the given boundaries are connected to z, or more 
precisely to a disk of radius e around z. Individually they are non-universal and furthermore go to zero as the lattice 
mesh size (or e) goes to zero. However, we find numerically and prove via conformal field theory that the ratio 

= P ™W , (1) 

y/P L (z)P R (z)U h 

is a universal function of z in the limit that s goes to zero, depending only upon the boundary conditions on the sides 
of the rectangle. 

For most of this paper (except the last section) the boundary conditions (b.c) are assumed to be open, or free, on 
the top and bottom sides, and either open or wired on the left and right sides. The behavior of C(z) near the left 
and right sides depends strongly on our choice of b.c. However, for rectangles with width W greater than a few times 
their height H, C{z) goes exponentially to a universal constant, Co, for points that are on the order of W away from 
the left and right sides, regardless of the boundary conditions on those sides. We find that the asymptotic value Co 
is the same as that found for the case of point anchors 0, 0] : 

o7/2 5/2 

C° = 3 3/4 r(1/3)9 / 2 - 1-0299268.... (2) 

This agreement is expected because for points z far from either vertical side the difference between anchoring to 
boundary points or small intervals becomes negligible. 

Furthermore, we find the surprising result that with wired boundary conditions on the left- and right-hand sides 
C(z) depends only upon the horizontal coordinate x and not on y, even though the individual functions Pl(z) etc. 
have a strong dependence on y. With wired boundary conditions all clusters touching the boundary are assumed to 
be connected together, so that if there is a crossing cluster all other clusters touching either boundary are also part 
of it. For this boundary condition, C(z) goes to 1 as z approaches the left or right sides, and remains within a few 
percent of 1 for all z, so that factorization is a good approximation everywhere in the rectangle. 

The results simplify particularly nicely for the case of a long rectangular system which we approximate as a semi- 
infinite strip (of unit width). In that case, we find 

nl , r ^(-1/2, -1/3, 7/6, e- 2 ™) 
C(x) = C - 



a / 2 Fi(-1/2,-2/3, 5/6, e- 2 ™) 
^o(l-|e-- + 5 ||e-- + ...) (for^oo) (3) 

where x is the distance from one end. We also provide numerical confirmation of this result. 

If we consider open boundary conditions on the left- and right-hand sides, then C(z) remains universal and ap- 
proaches Co far from the vertical sides. However, near the left- and right-hand sides, C(z) depends upon both x and 
y. Furthermore, when z approaches the left or right side C(z) goes to zero, so the factorization breaks down. 

We have also simulated C(z) with periodic boundary conditions on the horizontal sides. Here, of course, all P's 
are trivially independent of y and therefore C(z) is again a function of x only (for either open or wired b.c. on the 
vertical sides). C(z) goes to a constant value different from Co, viz. C\ w 1.022 ... at points far from the ends of the 
cylinder. We do not have a theoretical prediction for this value. 

Section |TT] gives our numerical results which show approximate factorization and interesting y— dependence with 
wired b.c. , and section HTT1 presents our theoretical derivation of those results; first for the case of the semi-infinite 
strip (including open b.c. results and complete expressions for the densities) and then for the complete rectangle. 
We compare these predictions with further numerical results. In section IIVI we consider the problem of finite-size 
corrections around an anchoring point. We show that our formulas for wired b.c. predict these corrections. In section 
|V]we present our numerical results for periodic b.c, and section fVTl gives our conclusions. The appendix presents a 
full derivation of the semi-infinite strip densities, which confirms that our expressions for Pl(z), Pr(z), and Plr(z) 
are exact, as is the y— independence of C{z) with wired b.c. . 



II. NUMERICAL RESULTS FOR C(Z) 



To investigate the probabilities Pl(z), etc., we carried out simulations in rectangular systems of dimensions 63 x 127, 
and 127 x 255. (The arrays used in the computer code were actually exact powers of two, but one column and row 



FIG. 1: Plots of Pl(z) (top), Plr(z) (center), and C(z) (bottom) for a 127 x 255 system with wired b.c. on the left and 
right-hand sides. 



were left open, along with periodic b.c., to efficiently simulate the open boundaries.) We used a square lattice and 
considered both bond and site percolation at the critical thresholds 1/2 and 0.5927460 @, H, respectively. The 
random-number generator used was R(471, 1586, 6988, 9689) given in [lfj. We kept track of the average density of 
clusters that touched the left, right, and both sides of the rectangular systems, where the density at a point is simply 
the number of times a cluster touching the desired boundary or boundaries includes that point, divided by the total 
number of trials. In Fig. Q] we show plots of the densities for bond percolation on a lattice of size 127 x 255 sites 
with wired boundary conditions on the left and right-hand sides. The top figure shows Pl(z); the plot of Pr(z) is 
identical but flipped horizontally. Along the left and right-hand boundaries Pl(z) is constant and equal to 1 and II/j 
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respectively, both a consequence of the wired b.c. In the intermediate region there is an exponential drop-off in the 
density. 

In the center plot in Fig. [1] we show Plr(z), which is roughly independent of x away from the ends. The lower figure 
shows C(z) defined by Eq. ([T]), and here one can see the striking result that C(z) depends upon the x-coordinate but 
appears to be independent of y, in spite of the strong y-dependence of the functions that define it. Note that the 
range of the vertical scale now goes from 1 to 1.03. 

At the two wired boundaries x = and x = w, C'(z) goes to 1. For bond percolation, where all sites are effectively 
occupied while the bonds are diluted (occupied with probability p), C(z) is identically 1 at these two boundaries 
because Pl(z) — > 1, Plr(z) — > IT^, and Pr{z) — > IT^ as x — > for a given y (and similarly for x — > w). For site 
percolation, where sites including those in the first and last columns are occupied with probability p = 0.592746 . . ., 
we have Pl{z) — » p, Plr( z ) —> pHh> an d Pr(z) — * pHh as i -> for a given y, so here C(z) — ► I but only on 
the average, not identically as in the bond case, and there are small fluctuations. We do not show the plots for site 
percolation as they are quite similar to those for bond percolation. 

Away from the left- and right-hand boundaries, C(z) approaches the value Co given in Eq. ((2]). The quantity Co first 
appeared the context of the densities of clusters touching one or two boundary anchor points, where after just several 
lattice spacings away from the anchors the analogous C(z) was found to go to Co everywhere Q. The reason that 
the same constant appears here is that, from a distance, the interval looks like a point. Furthermore, as discussed in 
more detail below, a conformal transformation converts the interval problem to the point anchor problem, and shows 
that indeed C(z) far from the vertical boundaries of the rectangle asymptotes to C . 



Fig.[2]shows a plot of the ratio Pl(z) / Plr{z) for the 127 x 255 system as shown in Fig. [TJ showing quite clearly that 
this ratio is only a function of x and not y for wired b.c. Likewise, Pl (z) / Pr(z) only depends on y (this follows because 
C(z) is only a function of y). Thus, all three quantities Pl(z), Pr(z), and Plr(z) have a common y dependence for 
a given value of x. We will use this to motivate our theoretical results. 

In Fig. [3] we show similar data for a 63 x 127 bond percolation system, this time plotted as contour plots. The 
apparent y independence of C(z) is quite clear. For much smaller system sizes (e.g. 15 x 31), finite-size effects do lead 
to visible curvature in the contours of C(z) (not shown here). 

The data for x — 1, . . . , 5 for a system of size 127 x 255 are shown in the upper plot in Fig. [U showing the fluctuation 
in the data about the constant values. These rather large fluctuations are due to the fact that the probabilities Pr{z) 
and Plr(z) are quite small for z near the sides. Here 1.14 • 10 9 samples were generated, a similar number as for 
the other plots. The lower plot shows the same simulation with open b.c. for comparison; the dependence on b.c. is 
readily apparent. 

In Fig. [5] we show C opcn (z) for the case of open boundary conditions on the left and right-hand sides. Here the 
behavior differs markedly. While C opon (z) still approaches the value Co away from the sides, closer to them a significant 
y-dependence of C(z) can be seen. Also, C(z) does not go to 1 at those sides as it does in the wired case, but rather 
drops to 0, highlighting the sensitivity of the factorization to boundary conditions. 
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FIG. 2: Plot of Pl(z)/P lr (z) for the wired 127 x 255 system of Fig.[T] 
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FIG. 3: Contour plots of Pl(z) (top), Plr(z) (center), and C(z) (bottom) for a 63 x 127 system with wired b.c. on the left 
and right-hand sides (similar to Fig. [T]but here as contour plots). In the plot of C(z) (lower figure), the contours are at 1.005, 
1.010, 1.015, 1.020, and 1.025. The lower figure may be compared with Fig. 1101 which shows a similar plot for open b.c. where 
the behavior is quite different. 



III. THEORETICAL RESULTS 

In this section we use boundary conformal field theory to calculate C(z). This quantity (see |T|)) is given by a ratio 
of various densities. Our work in @ illustrates the method by which such densities can be derived from correlation 
functions in conformal field theory with c = 0. The three densities we are interested in arise as different conformal 
blocks for the same set of operators, being distinguished by fusion rules, as we will see. Thus in the upper half-plane 
H, with w = u + iv, 

Plr(w), Pl(w), Pr(w) ~ (c6i,2(wi)01,2(u2)^(w 3 , Wz)(t>l,2{u4)(t>l,2{u 5 )) , (4) 

where the boundary operator 4>\^{u), with conformal dimension ftx,2 = 0, implements a change from open to wired 
(free to fixed) boundary conditions at the point u, and the "magnetization" operator ip(w) := 4>3/2,3/2(w) 1 with 
conformal dimension = h^p = 5/96, measures the density of clusters at a point w in the upper half-plane. The 
other quantity of interest, IT/j, is simply the probability that there exists a crossing cluster in the rectangle, given by 
Cardy's formula (which itself follows from a correlation function like (|4|) but without the operator) [l[ . 

These correlation functions can be constructed from the conformal blocks allowed by (U). However, computing these 
conformal blocks is a formidable task involving a six-point correlation function that depends on three independent 
cross-ratios. We can make the calculation tractable by exploiting the observed independence on the vertical coordinate 
in the rectangle: C(x, y) = C{x). Letting y 3 — > (so v 3 — > in H), the bulk magnetization operator ip in H is replaced 
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FIG. 4: Top figure: C(z) as a function of Y for X — 1, 2, 3, 4, 5 (bottom to top) for a system of size heightxwidth = 127 x 255, 
with wired b.c. on the sides. Bottom figure: same simulation with open b.c. (Note the change of vertical scale.) 



with a boundary magnetization operator, ^1,3 (1*3) so that 

Plr{u),Pl{u),Pr{u) ~ (</>l,2(ui)^l,2(«2)<Al,3(W3)^l,2(W4)^l,2(W5)) • (5) 

Note that ([5]) is a five-point function involving only boundary operators that we will compute shortly. In this 
section, we assume y-independence throughout. However, the appendix shows that our results are valid without this 
assumption in the case of the semi-infinite strip. 

We next determine ([5]) in the upper half-plane, and then transform the result into the appropriate geometry. (This 
transformation leaves C invariant [5[.) 

In subsections IIII Al and IIII Cl we find closed expressions for Plr(x), Pl(x), and Pr(x) in the semi-infinite strip and 
rectangle. The appendix verifies that for the case of a semi-infinite strip these expressions are valid without assuming 
y-independence. 



FIG. 5: Plot of numerical results for C op en(z) for a system of size 127 x 255 with open boundaries on all sides. In contrast 
to the wired case, Fig 

[2 Cop en ( z) goes to a value much less than 1 at the ends, going to zero as the system size increases to 
infinity. The dependence of C op cn(z) upon the y coordinate is evident here. 



Note that the variables z and x represent bulk and boundary densities in the rectangle or strip while w and u 
represent bulk and boundary densities in H. 

A. Semi-infinite strip 

We next consider the semi-infinite strip, which approximates the behavior near the end of a long rectangle. This is 
equivalent to replacing the right-hand interval with a point operator because it is vanishingly likely that two distinct 
clusters emanate from the distant end. Depending on which of the densities we are calculating we make use of the 
fusion rules to replace </>i,2(w4)</>i,2(m5) with either l(oo) or <j)i. 3 (po), which significantly simplifies the calculation. 

Because <j)i. 3 is the boundary magnetization operator, the fusion ^1,2 x 4>\.2 — $1,3 conditions the wired interval 
to connect to other objects marked in the correlation function. For example, to calculate Pr(u) and Plr(u) and 
then transform the result into the semi- infinite strip we insert the ^1,3(00) operator, because these two quantities 
require that a cluster connect the point at infinity with u. (Here and in the remainder of this section we ignore the 
scaling factor and operator product expansion (OPE) coefficient in this fusion since they cancel in the ratios that we 
calculate.) 

The fusion 0x,2 x ft 1,2 — 1 places no condition on the connectivity of the wired interval. So we get configurations 
with the interval isolated and with it connected to other points. Note that the identity fusion encompasses all of the 
configurations in the <p±^ bulk-boundary fusion as well. This mirrors Cardy's derivation of the horizontal crossing 
probability, where 11^ corresponds to the 4>i.3 conformal block while the identity block counts all configurations 
including those with a horizontal crossing [l|. 

As mentioned, one can calculate Pr(u) and Plr(u) using the correlation function 

(01,2M^l,2M^l,3M01,3M} = (U4 - U 3 )- 2 / 3 F ( (^- U ^~U 3 ) \ 

\{u 3 - Ul)(U 4 - u 2 ) J 

The null-state for 4>i^ gives 

= ( -, V3 + , V2 ~ ^ ^~ - |C) (01,2(«l)01, a (« 2 )^,3(«3)01,3(«4)) • (7) 

Applying ([7]) to © and using conformal symmetry to take {u\, U3, U4} {0, A, 1, 00} in the usual manner, so that 
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A is the cross-ratio, we find 



= F"(X) + 2 ^ '(A) 



r F(A) 



3A(1-A) v ' 9(1 — A) 2 
The solutions yield the conformal blocks 

Fi(A) = (l-A) 2 / 3 2 Fi(l/3,4/3,2/3,A), and 
F 013 (A) = A 1 / 3 (l-A) 2 / 3 2 F 1 (2/3,5/3,4/3,A). 

Here the index on the blocks represents the corresponding fusion channel of the ^ 2 operators. 



(a) 



(8) 



(9) 
(10) 



FIG. 6: Schematic representation of the two types of cluster configurations consistent with {<f>i : 2(0)^1, 2(A)0i, 3(1)^1,3(00)), 
translated into the strip geometry. 



z=x + iy 




K 
FIG. 7: Semi-infinite strip mapping. 



By the above, we may identify the configurations with (0, A) connected to 1 and infinity with the conformal block 
i<0 13 (A). Upon mapping into the strip this correlation will become Plr(x) as illustrated in Fig. [6ja). 

The conformal block -Fi(A) represents the configurations with (0, A) isolated from, as well as connected to, 1 and 
infinity. Thus upon mapping to the strip (see Fig. [7]) this block becomes Pr(x), which includes all configurations with 
x connected to R, as illustrated in Fig. [6^a&&). Note that i*i(A) (in H) includes configurations where u and R are 
not directly connected, but instead connect through the wired boundary conditions on the left side. This subtlety 
will be discussed in more detail in section UlI Bl 

Having found Plr{u) and Pr{u) we need only determine Pl{u). Configurations of this type don't need to cross to 
the right-hand side, so instead of placing 01,3 at infinity we place an identity; thus Pl{u) is proportional to a three 
point function. 

Thus, including the proper normalizations and scaling factors we find 



Plr(u) 
Pl{u) 

n,, 

Pr{u) 



h,2 (0)^1,2 (A) [l i3 ] 01,3 (1)01,3 (£)) 
^2(O)01, 2 (A)0i, 3 (l)) 



Mo)0i, 2 (A)0i, 3 (£)> 

h,2 (O)01, 2 (A) 01,3 (1) 



C112C222C 2 ^ 3 -p0 1:3 (A) , 
CmA^fl -AJ-V3, 

Cii 2 A 1/3 r 2/3 and, 

r 2/3 ^(A) , 



(11) 

(12) 
(13) 
(14) 



where the variable £(— > 00) indicates the order of magnitude of the length of the rectangle, the bracketed subs crip t 
notation denotes the fusion channel propagating between the operators separated (following Bauer and Bernard [111]), 
and the indices on the (boundary) operator product expansion coefficients djk follow the convention of jl2j with 
index j representing the j— leg boundary operator 0i,i + j. Both operator product expansion coefficients that appear 
above are previously known: C112 2 = 27r\/3/r(l/3) 3 = 0.56604668 ... is most familiar as the leading coefficient in 
Cardy's formula []0] , while C222 = C'o arose in the factorization formula in [j| @ . Note that the expressions for Pl (u) 
and could also be obtained as limiting cases of Cardy's formula when one interval shrinks to a point. 
For future reference, note that if the interval L is not connected to cither u or R we have 



Plr{u) 



"1,3 



(l)01,2(A) a[ l,2]01,2(Q)01 ) 3(O)=^ 1 , 2 r 2/3 Ga 01 , 2 (A) , 



(15) 
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where 

Ga0 1 , 2 (A) = (l-A) 2 / 3 2 F 1 (l/3,4/3,2,l-A) and = gpl^ • (16) 

This can be deduced by observing that Plr(u) = Pr(u) — Plr{u) and using the crossing symmetry relation for this 
correlation function, 

^30i, 2 G,90 li2 (A) = F 1 (X) - Cn2C222-F 1 1 , 3 (A) • (17) 

In (|15[) we label the propagation channel and conformal block with d4>i,2 m order to be consistent with the underlying 
logarithmic conformal field theory However, for present purposes the fusion 0x,3 x 0i,2 is more complicated to 
interpret physically than the 0i ; 2 x 0i,2 fusion discussed above, so while we include the notation for accuracy, the 
distinction is unimportant here and in what follows. 
Combining fT] and ([^-([H]) we find 

C(X) = CoVT-x ^(W/^aa) . (18) 

^1(1/3, 4/3; 2/3; A) 
Now map from the semi-infinite strip S = {z = x + iy\0 < x, < y < 1} via 

cosh 2 (7rz/2) 

W M = T27 7^7 » ( 19 ) 

cosh (ttx/2) 

which takes the point x on the lower side of the strip to 1 in the upper half-plane (see Fig.[7|). 
The cross-ratio is given by 

A = w(0) = sech 2 (7nz;/2) . (20) 
Using (|18[) and standard hypergeometric identities involving quadratic arguments [l3[ gives 

C( X ) = C fi(-V2,-l/3,7/6,e-^) 
^(-1/2,-2/3, 5/6, e~ 2 -) 

The appendix shows that (|2ip is in fact exact, and follows without the assumption of y-independence. 

Now far from the finite end of the strip C(x) must reduce to the analogous result involving the density of percolation 
clusters constrained to touch two isolated boundary points [5j. Thus Co = C222 and ((2|) follow. Alternately, if x = 
then (as mentioned) the wired boundary conditions ensure that Pl{z) — 1 and Plr(z) = Tlh — Pr.{ z ) so that 
C(0) = 1. It is a non-trivial check of (|2"Tj) that these limiting cases are indeed consistently recovered. 

Expanding (f2"Tj) . we find 

C W _«,(!- | e - + Jie- - ^e-- + ) , ) - (22) 

In Fig. [5] (upper curve), we compare this prediction with simulation results from a system of size 63 x 255. Fitting 
the numerical results to a polynomial in s = exp(— 2nx), we find C(x) — 1.02953 — 0.059993s + 0.033494s 2 . . . = 
1.02953(1 - 0.058272s -I- 0.032933s 2 . . .) while (J22J) gives C (l - 0.0571429s + 0.0333267s 2 . . .), so the agreement is 
quite good. 



B. Density in the semi-infinite strip 



In the previous section, based on numerical evidence, we assumed that the ratio C{z) in the semi-infinite strip 
was independent of vertical position and performed the calculation in the simpler case when z is on the boundary. 
However, the simulations further suggest that any ratio of two of Plr(z), Pl(z), and Pr(z) is also independent of 
y, as shown in Fig|2j Given this stronger condition (which, for the semi- infinite strip, is proven in the appendix), 
knowing any one of these three functions immediately determines the other two via the expressions in section IIII Al 
In this section we exploit this idea to find Plr(z), Pl(z), and Pr(z) and also the corresponding quantities for open 
b.c. 
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FIG. 8: C(x) vs. s = e~ 2lrx for wired b.c. on the vertical sides, for open b.c. (upper curve) and periodic (lower curve) b.c. on 
the horizontal sides, for systems of size 63 x 255 (open) and 64 x 128 (periodic). Here x = (X — l)/63 (open) and (X — l)/64 
(periodic), where X is the lattice coordinate of the column. Data are fit to fifth-order polynomials as shown. For the open case, 
the first few terms compare favorably with the prediction, Eq. (|22[) . For the periodic case, we have no theoretical prediction. 



In the semi-infinite strip we can easily calculate Pl (z) based on results for the density of clusters anchored to an 
interval in [5j. As the side R is taken to oo, the probability of spanning the length of the strip is negligible compared 
to the probability of connecting to side L. Thus there is no distinction between wired and open boundary conditions 
on R, and we can equate Pl{z), with p^(z), the density of clusters attached to the left side when the sides have open 
boundary conditions. (In what follows we differentiate between the boundary conditions for our densities by using £ p' 
to refer to open sides and 'P' to refer to wired sides.) 

The density of clusters at a point w = u + iv in the upper half-plane that are attached to an interval / = (u\, 112) 
is 1 

\2vl - V J {w - ui)(w - U2) 

We next transform this into the semi-infinite strip using the mapping (|19p . The points {i,0,z} in the strip map to 
the points {0, A, A cosh 2 (f (u + iv)) } in the upper half-plane so that the corresponding cross-ratio is 



_ (A-0) (Acosh 2 (f(x-iy)) - Acosh 2 (|(x + iy))) _ /sinh(Trx) + i shuVy) \ 2 
''~ (Acosh 2 (f (ir-Hy)) -0) (A cosh 2 (f (a; - iy)) - A) ~ ~ \sinh(7ra;) - i smfry) ) 



This leads to 

11/ 



P l{z )= Pl{z) -L ( 2 Sin ^ _ ) (25) 

(sinh (nx) + sin (Try)), 



The expression for Pl(z) then determines the other two densities. Using (fTTj) . ([12]) . and (jT4j) and including all 



11 



relevant transformation factors we find, for points x on the side of the strip, 

1 /3 

P LR (x) = U h C 222 (jz^) tFx (-i.-i.Z.e- 9 ™) . (26) 

T^=^) > and ( 2? ) 



^ = fl£ W^l-e-^) ) ^ 1-3-2' 6' 6 J < (28) 
where A = sech 2 (7ra:/2), u/(a;) = 7r\/l — A, and hypergeometric identities have been used. Thus 

n h P L (z)C 222 e™/ 3 / 1 17 _ 27ra 



n,P L (z)e 2 -/3 ^ 2 15 _ 27rx 



*w = -; 2 ' 42/3 ^^-3.-2' 6' 6 J < (30) 

where Pl(z), given in (|25|) . contains all the y dependence in these expressions, and il^ is given in (|13p . 

These densities apply with wired boundary conditions on the left side, but the analogous expressions can be found 
for open boundary conditions. One more result is required: pr(z), the density of clusters attached to the distant 
right-hand side when the left side has open boundary conditions. In the upper half plane this quantity is given by 

p R (w) = (tP(w, u#i, 3 (0> = v- 5/4S L - J ) , (31) 



2k -?| 2 



with £ — > 00, which becomes 



(32) 



U h I 2, ■ 1 2/ \ . • 2 / xx\5/96 



(4tt) 1 /3c 1 
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(^ 2 (sinh 2 (7rx) + sin 2 (^))) J/ (sinh(Trx) sin(7ry)) 11/48 (33) 



n ' lPi if) sinh^^vrx) (sinh 2 (^) + sin 2 (vry)) 1/6 , (34) 



c 112 2 41/3 



when transformed into the strip. The last equation is included to simplify comparison with (|30|) . The normalization 
of pr(x) is consistent with that chosen for (|25|) (the same non- universal factors are omitted). A check of this is that 
Pr(z) —pn{z) — > as x — > 00, which we expect since the boundary conditions on the left hand side are unimportant 
in that limit. 

With this result we can calculate the density of crossing clusters, Plr{z) with open boundary conditions on the 
left hand side. Note that the density of clusters that attach to the distant right side but not the left is independent 
of the particular boundary conditions we place on the left side, and as discussed the boundary conditions on the 
right-hand side are unimportant in the strip. This implies that Pzr(z) = Pzr(z), with Pzr(z) — Pr(z) — Plr(z) 
and Plr( z ) — Pr( z ) ~ Plr{ z ), where, e.g. Pir(z) is the density at z of clusters that connect to R but not to L. An 
explicit result for this quantity 

Plr(z) = P L r{z) = 1C n^/S^ 2 2*1 -'O' 3 ' 1 " 6 ( 35 ) 



16e- 2 ™/3 Cn2 2 l \3'2 

follows from (l29l) and (l30l). and leads to 



, , n h P L (z)sinh 1/3 ( TO ) Msinh 2 (^x) +sin 2 (7rj/)) 1/6 
PLr[z) = I 4V3 (36 ) 

i4T a012 tanh 5/3 (7re) ^ ^4 5 n 



after simplification using hypergeometric identities. 



2 Fi ( -,~,2,tanr/(7ra) 
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Thus for a strip with open boundary conditions we have 

Plr(z) 



Copenf^} y) — 



\JPL(z)p B {z)T\ h 
sinh 1 ^ 6 (7ra:) (sinh 2 (7ra;) 



snr 



1/12 



ifa0 12 sinh 1/6 (7ra;) tanh 5/3 (7ra;) 



di2 2 5 / 3 (sinh 2 (7rx) + sin 2 (Try)) 



1/12 



1*1 



, 2, tanh 2 (7rx~ 



(37) 



This expression is considerably more complicated than (|21[) . and exhibits a dependence on y that is not present with 
wired ends, as shown in Fig. [5J In addition < C open (x,y) < C222 while 1 < C(x) < C222 so that factorization is a 
much less accurate approximation near an open anchoring interval. 




FIG. 9: Theoretical prediction for C(z) (upper, red; Eq. (|21[) ) and Copcn(z) (lower, blue; Eq. ((3TJ) ) near the end of an infinite 
strip, which are comparable to the simulation results shown in Figs. [1] [4] and [5] The two plots differ only in scale: the right-hand 
plot focuses on the range of C(z) which is obscured on the left. 




FIG. 10: Simulation results (above) of C op en(z) on a rectangle of size 63 x 127, compared with the theoretical prediction of 
P?)) for a strip. Contours are at C = 1.029, 1.025, 1.02, 1.01, 1.0, 0.95 and 0.9, going toward the ends. 



To compare this result with simulations, we considered a system of size 63 x 255 with open b.c. on all four sides, 
generating 10 8 samples. In Fig. 1101 we show the contours predicted by Eq. (|37l) (lower figure) and compare them with 
the contours of the numerical simulations of the 63 x 255 system (upper); clearly the behavior is similar. This may 
also be compared with Fig. [31 which shows the contour lines with wired b.c. 
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To make a more quantitative comparision, in Fig. [IT] we compare the results at one end of the rectangular system 
with the theoretical strip results, Eq. (|37[) . Asymptotically for large x, Eq. (|3T|) behaves as 

<W*. V) = Co (l - ^3 )6ll + 5 2To (2 ^ ) e ' 5W3 + 13 + 7 4 7 (2 ^ e— + 0(e-"»/»)) (38) 

The first correction term implies that the contours of constant C(x,y) are determined by x = xq + 3/(5%) log[ll + 
5 cos(27ry)] where xo is a constant. This is roughly sinusoidal and has an amplitude of 3/(5tt) log(8/3) w 0.1873, which 
is consistent with the behavior seen in Fig. [lOl Thus, the sinusoidal behavior in the y-direction of the correction term 
for the open b.c. persists for all x, although its amplitude drops off exponentially in the x-direction. 
Along the centerline y = 1/2, Eq. (j38|) yields 

c opcn (x,i/ 2 ) = Co (i - 4 ^ r 22 ( /C e " W3 + r~ 2nx + ( e " ll7rx/3 )) 

= C - 0.33492472 e ~ W3 + 0.14713240 e~ 27TX + ©(e^ 11 ™/ 3 ) (39) 

In Fig. [TlJ we plot ln(Co — C opcn ) vs. x, which should approach a straight line with slope — 5ir/3 for 0<Ci« w/2, 
where w — 4 is the width of this system. The agreement between the prediction and theory is seen to be excellent for 
x < 1. 

The first two exponents in |39|) arc close to each other, but far from the next exponent (— 117t.t/3). In Fig. [Tl] 
we also show how the numerical data compares with using those first two corrections terms in (I39p . They give an 
excellent fit for 1/2 < x < 1. 
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FIG. 11: Results of ln[Co — C op en{x, 1/2)] vs. x along the centerline y — 1/2 for a lattice of size 63 x 255 (data points) compared 
with theoretical prediction for the half- infinite strip, Eq. (|37p (red solid line), one term (blue dashed line) and two terms (green 
dotted line) in the asymptotic expansion, Eq. (f39|) . 



C. Rectangle 

We now consider the finite rectangle. Here, we determine the five point function in ([5]), and extend our results to 
y > by assuming the y-independence observed numerically. 

The Kac null states and associated differential equations are less helpful in the case of correlation functions with 
more than four operators. Instead we find the conformal blocks using vertex operators in a Coulomb gas formalism 
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using the methods described in [l^|. This method, as it relates to four-point function, is fully described in [lij. The 
generalization to five-point functions is straightforward. 

The central charge of the conformal field theory is c = so that the Coulomb-gas background charge is — 2ao = 
— 2/V24- The vertex operators with charge a, V a (z), correspond to conformal operators with weight h — a(a — 2ao). 
Thus for (f>i r 2 with hi t % = we require a = or 2«o; for ^3 with hi 3 = 1/3 we have a = — 2ao or 4ao- By the 
charge neutrality condition only correlation functions with net charge 2«o are non-zero. In the plane, these take on 
the form 

<n^(*)>=nte-*) 2o<oy - ( 4 °) 

• i<J 

In order to calculate correlation functions consistently we also need to include zero weight screening operators. These 
non-local operators are formed by taking the integral of weight 1 operators around a closed loop; they compensate for 
any excess charge in the original set of operators. Since h = 1 allows either a — — 4ao or 6«o we have two screening 
operators at our disposal. 

j> duV ( i ±5)ao (w) , (41) 

where T is the closed contour of integration, which must non-trivially wrap around the singular points of the correlator. 

For our purposes, in the half-plane coordinates, the simplest expression for ((SJ with the correct scaling weights and 
charge neutrality is 



(0i, 2 (-oo)0i, 2 (l - 1/M)<k2 (0)^1,2 (^)<k3 W) = f (^aao (-00)^200(1 - l/^V 2ao (0)V 2ao (iy)V- 2ao (l)V- iao (u))du , 

(42) 

with v, (j, € (0,1). The intervals (— 00, 1 — 1/a 1 ), and (0, i/) are wired, while the point at 1 is on the free boundary 
outside of these intervals. 

In most cases the contour T can be chosen so that it can be replaced with an open contour between singular points 
of the correlation function [l4j. This means we can rewrite the correlator as an integral 

Ki-m)(i-m + ^)V /3 [( «-i V /3 dW; (43) 



— ) J 7 \u(u — v)(u— 1 + 1/ /i) 

with 7 a segment of the real line between two of the five points. 

The various choices for 7 lead to the following formulae, where in each case we pick the branches of the integrand 
in order to obtain a real solution: 

IYl/Sl 2 v 1 ! 3 / \ 

(— ,i-i/rt- 7 *w« = r ( ^ {1 -, mi - ll + l ^ T '(i*-i--i^^) < 44 > 

(l-l/*0)= 7 WW- AM (i _ y)i/3(i _^ + ^ )i/3 Fl (i ;i> - §;i |_^ >1 - M ) (45) 

(0,^=7^^ = (1 ^)i/3 F i (I; §1 llrfe. T-fe) (46) 



27 

(l,oo) =7 =► = 1 =(f(1-/')(1^) 2 (1-/' + H) F x (|;| ! |;2| I /,1-m + ^) (48) 



27 

These expressions, which involve the Appell hypergeometric function Fi(a; bi, b%) c\zi, 22), are conformal blocks 
of $5§. The blocks represent the configurations in Fig. [12] as indicated by their subscripts, i.e., Wab^Bi represents 
AUB1UB2, while Wb x represents the B\ configurations only. Thus, for example, Wa — Wab x — Wbx — Wab 2 —Wb 2 - 

In order to justify our assignments for these conformal blocks we take the limits v or /1 — > for comparison with 
those found directly from the five-point function. Taking the leading term of the five-point function as either v or 
(j, — ► gives the three or four-point functions of section IIII Ai which allows straightforward identification with the 
physical configurations. 
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(A) 





FIG. 12: The three types of configurations resulting from wired boundary conditions in the rectangle. 



z=x + iy 




1-1/u. 



FIG. 13: Mapping the rectangle onto the upper-half plane: the left side maps onto (0, v), the right side maps onto (— oo, 1 — 1/fi) 
and the point on the bottom side maps to 1. 



As v — > the weights of the configurations in Fig. [12] are 

Oi -» Cii 2 I/ 1/3 (0i,2(-Oo)0i, 2 (l-l/At)[l,3]^l,3(O)^l,3(l)) =^ /3 C Y 112 2 C222^ 1 , 3 (m) 



(49) 

i,2(-oo)0i, 2 (l - l/^)0i,3(l)> = CW 1/3 , (50) 
Vb 2 CW 1/3 <^i,3(l)0i,2(l - l/A*)s[i,a]0i^(O)^i,2(-oo)) = ^ 1/3 Cii 2 ^ li2 G^ li2 ( M ) , (51) 



where the notation f2 is used to emphasize that these expressions come from the original correlation function and 
therefore may have a different overall normalization than our conformal blocks. These limits are easily deduced by 
inserting the leading order term in the boundary operator product expansion for 0i, 2 (O)0i j2 (i/) into ([5]) as discussed in 
section llll Al if the interval (0, v) (the left side in Fig. [T2"f is isolated from the other structures the leading contribution 
will be 1(0), and if it is attached to the other structures the leading fusion term will be Cn2^ 1 / 3 0i,3(O). 
For comparison we take the limit v — ► in the expressions (|44"]) - (f48|) for our conformal blocks. We find 



W AB2 
Wb 2 



r(l/3) 2 vV* 
r(2/3) (l-A^/3 



Fi 



'1. 2 _2, 2 

,3' 3' 3' 3 



0,^^ 1/3 fiM, and 



-in 



(Ml - fi) 2 ) 1 ^ 3 Fi (|; §, §; 2|0, 1 - /x) = -^=i// 3 G 901 » 



4tt 
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(52) 
(53) 



Taking into account the crossing symmetry (|17p we see that Wab 2 ~ + and Wb 2 ~ ^b 2 , validating our 
labels for these two blocks. We find that (|4"S")) . (f4"6")) . and (|4"T)) are C(^°). Thus these three blocks must represent the 
B\ configurations. In order to specify the inclusion/exclusion of A or _B 2 in these blocks we next perform a similar 
analysis for ji — ► 0. 

In order to take the limit fi — > we must use the expansion around negative infinity; lj2 (— £)0i, 2 (1 — 1/a 1 ) — ► 1(— £) 
or Cn 2 /i 1 / 3 ^ 2//3 (/)i i 3(— £) with £ ^> 1/fj, and eventually £ — > oo. Here 



^ - Cn2fi 1/3 e /3 (4iA-O0iAO)[iMiA^iAl)) =M 1/3 Cii 2 2 C 222 F 013 H , 
r> Bl - C 112 m 1/3 £ 2/3 <^i,3 (-6^1,2 (0)^,2 Ma[i,2] ^1,3 (1)) = M 1/3 Cii2#d%, 2 G a01 , 2 H , 

tl B2 -» (01,2(O)01, 2 M01,3(l)) = Ci 12 Z/ 1 / 3 (l-^- 1 / 3 . 

We then take the limit /Lt — > in the expressions (|i4] ) -(|48 |) . which yields 

\2 



w ABl 
w Bl 



r(i/3) 2 fi 1 / 3 

T(2/3) (1 - i/)V 

4?r / - x 2 \l/3 



I73 Fl (I 



• 2 _2. 2 

'3' 3' 3 



47T 



(^(1 - ^) 2 )" /J F x (|; |, |; 2|0, 1 ^) = ^L=^ 3 G d ^(v) 
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(54) 
(55) 
(56) 



(57) 
(58) 



The comparison validates our assignment of Wa Bi and W Bl . 
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Having identified these four blocks is sufficient for our purposes but for completeness we can also justify our label 
on the fifth block, WabiB 2 - If we deform the contour of integration from 1 — l//i to into the upper half plane so 
that it runs just above the real line from 1 — l//i to — oo and from oo to then we can express Wab^b 2 in terms of 
the other four conformal blocks, 

W ABlB2 = -e^/ 3 W AB2 - e 2m / 3 W ABl - e im / 3 W Bl - e 2 ™' 3 W B2 . (59) 

Taking the real and imaginary parts of this relation gives: Wab x + Wb 2 = Wab 2 + Wb x , a non-trivial check of our 
established labels, and Wab x b 2 = § (Wabi + W B2 + Wab 2 + which validates our final label. 

Now that we have the necessary results in the upper half-plane, we transform them into the rectangle. The analytic 
function that takes the rectangle z = {x + iy\x £ (0, R), y £ (0, 1)} to the upper half-plane mapping (i, 0, x, R, R + i) 
(with x £ (0, R)) onto (0, v, 1, — oo, 1 — l//x) can be written 

nc 2 (^(l-a)|q) 
nc z {xK{\ — a)\a) 

where nc(z|m) is a Jacobi elliptic function and K(m) the complete elliptic integral of the first kind. We determine 
the elliptic parameter a from R in a slightly non-standard way. The usual elliptic nome in this geometry would be 
q = e _7r / fl , but we prefer to expand in terms of q' = e _7ri? . Since R — > 1/R is equivalent to a — > 1 — a we find that 

using the standard expression for the inverse nome and elliptic theta function identities [HI . 
Using elliptic theta functions, the undetermined parameters in the mapping are 

ffn> 2, T.H M s ^(O^T^^Ag') 2 , 

v = /(0 = cn (xK(l - a)\a) = , and 62) 

1/4(0, q'Yv2{n^x/2, q'y 

= 1 1 MQ, gQ^MjWg, q 'f 

M l-f(R + i) l + j^cn 2 (xK(l-a)\a) d^{Q,q') 2 d^x/2,q'y 1 ' 
Thus the probabilities on the rectangle are 

Plr{x) = f'{x) l ^W A , (64) 

P L (x) = /'(x^Wab, , (65) 

Pr{x) = f{x)V 3 W ABl , (66) 

where the conformal transformation introduces the factor 

a/3_/V n 80(^(1 -q)| a ) V /3 _ / . K ( M^Y \ ^H^MWV) ^ V3 



"(j 2 ^ 1 a) nd(^(l-a)|a)J " (, * K {^4)*) Mi™MMi™li<t)J ' ^ 

We include this for completeness even though it cancels in the ratios that follow. 

We emphasize that these probabilities are not physically normalized, since certain non-universal constant factors 
related to the point operator in the correlation function are ignored. Note also that in arriving at the above expressions 
we have used the Jacobi elliptic functions cn(z|m), sc(z|m) and nd(z|m), and subsequently replaced them with their 
elliptic theta function equivalents. 



fi *- Mi** Ml 



- 1 • .. : . 



150 



FIG. 14: Ratio of measurements of C(z) (averaged over all y) to theory (|75[) for bond percolation systems of size 63 x 127 (left) 
and 127 x 255 (right) with wired b.c., plotted as a function of the lattice coordinate X. 



The final expressions for the probabilities of the various configurations are 



Paub 2 = 



p 



AUB1UB2 



^3(0)^! 02 1? 4 (^) 

^^ 2 (O) 2 04(O) 2 ^3( 1 f I ) 3 

r(f) U 3 (o) 4 ^i(^)^2(^)0 4 (^) 



1/3 



12 22 

3' 3' _ 3' 3 



2 (O) 4 2 (O) 2 tf 2 (W 



Elf) 

:§) 



1/3 



12 22 

3' 3' _ 3' 3 



*3(°r 3 (o) 2 3 W 
MO) 4 Mo^i^Y 



PauB! = 



47T 
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/ ^^(0)^ 4 (0)2^ 2 (i££) 3 

' l v 3 (O)%l(^)0 3 ( i f ; )04( i f i ) 

/ ^^ 2 (O) 2 04(O) 2 ^ 2 (T) 3 

^ 2 (O) 4 3 (O) 2 ^i (^)°^ 4 (T) 
l1 9 4 (O) 6 ^ 2 (^) 5 3 ( 



1/3 



12 2 2 

3' 3'~3' 3 



^3(0) 4 ' ^3(0)2^3 (1^) 2 J 



^3(0) 4 ' tf 3(0 )2tf 2 



17TX 1 
2 / 



1/3 



17TX 1 

4 I"/ 



■ ; V 2 

/ ■ 1 K \ l/3 



12 2 2 

3' 3' _ 3' 3 



2 (O) 4 2 (O) 2 ^ 3 (^)' 



/ 77. r \ 
2 ) 



F [LI 2 -2 

1 3'3'3' 



^s(O) 4 ' tf 3 (0) 2 tf 2 (HS) 



2 (O) 2 ^i (^) 2 tf 3 (0) 2 tfi 



(68) 
(69) 
(70) 
(71) 



^ 4 (O)203 tf 4 (0) 2 tf 2 



, 2 / 



47T 



^2(0)^3(0)^1 (^) 5 M^) 5 y /3 F (i.n. 2 

27\ *tf 4 (0) 6 i? 2 (^) 5 tf 3 (^) 5 / 1 V 3 ' 3 ' 3 ' 



^ 2 (0) 2 i? 4 (^) 2 tf 3 (0) 2 tf 4 (^)^ 



^4(0) 2 1?2 (T) 2 ' ^4(0) 2 l?3 (T) 2 , 



(72) 
(73) 



Combining (f68j) . ([7T|) . (|72]) and including the horizontal crossing probability 

27T3V3 / ^ 2 (0) 4 \ 1/3 A 2 4 z? 2 (0) 4 
r(±) 3 W0) 4 J 2 1 U'3'3'0 3 (O) 4 , 

given by Cardy's formula with cross-ratio [iv j(\ — fi + /^^), we find the main result of this section 

C(z) = C(x) 



-PauSi — Pbi 



This result is compared with simulations in the next section. 

Note that our results give exact expressions for P^(x), Pr(x), and Plr(x) in an arbitrary rectangle. 



(74) 



(75) 



D. Comparison with numerical results for a rectangle 



In Fig. [T3] we show the ratio of the measured values of C(x) to the theory Eq. (|75|) . for systems of size 63 x 127 
and 127 x 255. There is a slight overshoot near the boundaries, which was lowered to some extent by assuming that 
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FIG. 15: ln[Co — C{x)\ vs. x, wired b.c, for systems of size 63 x 127 (top) and 127 x 255 (center), and theory for a system of 
aspect ratio 2 (bottom). 



1.1 




v 



FIG. 16: Numerical values of Pl(x,d)/ Pl(x, 1/2) as a function of vertical coordinate y for X = 3, 4, 8, 16, 32 (top to bottom) on 
a rectangular lattice of size 31 x 63, with wired vertical boundaries (data points). We set x = (X — l)/62 and y = (Y — l)/30. 
The corresponding results for Pr and Plr are identical, as are the results for x replaced by 1 — x. The curves are drawn 
assuming the strip results (|25|l . which are are seen to describe the rectangular results quite closely. 



the effective boundaries of the system were 0.2 lattice spacings inside the system. There is a small undervalue in the 
center which decreases as the system size increases. 

Another way to illustrate the difference is shown in Fig. [15l in which we plot In (Co — C(x)) vs. x for a systems 
of aspect ratio 2 for two different sizes. It can be seen that as the size increases, the deviations from theory around 
x = 1 decrease. Thus, overall, we find very good agreement between Eq. (|75p and simulations. 

In Fig. [TB] we explore the question of the vertical dependence of the correlation functions. In section [TT] we observed 
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FIG. 17: Theoretical density plots of C(z) around a semicircular anchor of radius e = 1 in a half-plane from Eq. (|76[) . Contours 
are at 1.029, 1.028, 1.027. . . going inward. 



FIG. 18: Measured density contours of C(z) on a lattice system of size 127 with anchors at the top and bottom, for bond 
percolation. Contours around the anchors (going inward) are at 1.029, 1.028, 1.027. . .; scattered contours near left and right 
sides are at C — 1.03. 



that the various functions Pl(z), Pr(z), Plr(z) all have the same y-dependence, given by Eq. (|25|) for the case of a 
strip but unknown for a rectangle. In Fig. [THlwe plot the measured values of Pl(x,i/)/ Pl(x, 1/2) as a function of y 
for fixed values of x. Here we are normalizing by the value half-way up, P(x, 1/2); we find that the strip result ([25]) 
is indeed quite accurate here. The results for the other functions Pr(z), Plr(z) are identical. 

Note that the system displays an interesting horizontal symmetry. The y dependence of, for example, Pl{z) shown 
in the top figure of Fig. Q] at x and W — x is the same. Thus since Pn(x,y) = Pl(W — x,y) this implies that 
Pl{z) I Pr(z) is only a function of discussed above. 



IV. FINITE SIZE CORRECTIONS AROUND AN ANCHOR POINT 



In this section, we use some of the results obtained above to explain finite-size corrections observed in a different but 
related problem. In [5| we examined factorization of correlations in rectangles (or the half-plane) that are anchored 
at points on the boundary. The result is closely analogous to (|TJ) , except that the clusters are conditioned to touch 
specified points on the boundary rather than intervals. Here, everywhere except near the anchor points, C takes on, 
as mentioned, the same constant value Co found above. Further, C = 1 at the anchor points only, and most of the 
deviation of C from its asymptotic value occurs within of the order of 10 lattice spacings of the anchors. Since the 
deviations are only significant for a finite number of lattice spacings, they are a finite-size effect and vanish in the 
field theory limit. 
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FIG. 19: Comparison of (|76[) with simulations for C(r,0) along three radial directions 8 away from an anchor point, where r 
is the radial distance, and e = 1.34 lattice spacings. 



We can nonetheless use our conformal field theory results to model these corrections to C near the anchor points. 
We model the boundary anchor point as a semi-circle of radius e centered on the origin in the upper-half plane and 
place the other anchoring point at infinity. This mimics the case where the anchoring point has an effective finite size 
that is much smaller than its distance from other significant points and boundary features. The mapping w{z) — ee* 2 
from the semi- infinite strip to the model geometry implies the relation e/r — er' KX between the horizontal strip 
coordinate, x, and the radial coordinate r 6 {e, oo} of the target space. Inserting this relation into (|2l"j) we find 

C(r) = Co [ 2 ' , (76) 



P (_1 _2 5 e±\ 
2 ^!V 2' 3' 6' r 2 > 

with C independent of 9 (the contours are semi-circles). For large r, one has (similar to (|22p) 

C(r)=Co ( v 1 -^^ + 2^025^- 734825^--j " (77) 

so that the finite-size corrections drop off to first order as 1 /r 2 . This prediction is illustrated in Fig. [T7] 

To test this prediction, we measured C(z) for a system of 127 x 127 lattice spacings with the anchors at the centers 
of opposite sides, and free boundary conditions elsewhere, for critical bond percolation. The contours are shown in 
Fig. [TH1 In Fig. [TH] we compare ln[Co — C(z)] as a function of distance from the anchor point with (|76[) in the three 
directions indicated. The numerical data is scaled so that e = 1.34 lattice spacings in order to match the theoretical 
predication. There is an additional scale factor related to where to put the first point (the boundary location); that 
is chosen at 1.15 radii. These two scale factors have considerable leeway - for example, choosing 1.2 for e gives a very 
good fit for the perpendicular C(r), but not as good for the parallel (side) C(r). There is more noise along the side 
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FIG. 20: Plot of ln[Ci — C pC r(x)] for a periodic system of size 64 x 127, bond percolation, with Ci = 1.022, for X = 1, ... , 50. 
The equation for a linear fit is shown in the figure. 



(x — 0), because the probabilities are lower there. Also, as r increases, Co — C(r) decreases, so the noise grows. For 
the 45° case the two scale factors are multiplied by \[2 to be consistent. 

Note that, by the conformal transformation above, lnr in Fig. 1191 corresponds to x in Fig. 1151 for small x, since for 
small x the rectangle may be approximated by a strip. 

Thus, the overall agreement with theory is good, and the extent of measurable finite-size effects (of the order of 10 
lattice spacings) agrees with our previous observations. 



V. PERIODIC SYSTEM 



Next we consider a system with periodic b.c. on the horizontal sides and wired b.c. on the vertical sides. In this 
case we have no theoretical predictions for the factorization but just present the simulation results, for a system of 
size 64 x 127. We find that away from the boundary, C pel (x) approaches a constant C\ ~ 1.022. Further, C per (x) 
does not deviate from C\ by more than 2.2%, so factorization is a good approximation everywhere in the system. 

We determined the value of C\ by plotting ln(Ci — C por ) vs. x and adjusting C\ to get a good linear fit in the 
intermediate regime between the vertical walls and the center at x — 64. This plot is shown in Fig. [201 The slope of 
that fit w —0.09557 is reasonably close to the value — 27r/64 = —0.09817477 one would get from the expected leading 
behavior exp(— 27r(A — l)/64). Assuming higher-order terms are powers of this exponential term, in analogy with 
the case of wired b.c, we plot Cp er (a;) vs. s = exp(— 2tt(X — l)/64) in Fig. [5] (lower curve). We get an excellent fit 
to a polynomial with the result C pe r(a;) = 1.02200(1 - 0.042758s + 0.014994s 2 . . .). Thus, while the behavior of the 
periodic system is quantitatively different from that of the open b.c. case, qualitatively, the two are quite similar. 
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VI. CONCLUSIONS 



In this paper, we consider two-dimensional systems at the percolation point. We investigate correlation functions 
between a point in the bulk of a rectangular system and one or both of the two vertical boundaries. The b.c. on the 
horizontal sides are assumed to be open (except in the last section where we consider periodic b.c), and the vertical 
sides have either wired (fixed) or open b.c. We discover numerically that for wired b.c. on the vertical sides, all the 
correlation functions have an identical y-dependence, which implies that the ratio C(z) (see (fT[)) is independent of y. 
Further, this ratio is close to 1 for all z in the rectangle, so that factorization of the correlation functions is almost 
exact everywhere. 

Assuming y-indepcndcncc, we use boundary operators in conformal field theory to find C(z) and the correlations 
that comprise it everywhere in the semi-infinite strip (long rectangle) from the much simpler calculation of the density 
at the wall. Explicitly, we solve the differential equations that arise via conformal field theory when the bulk point 
has been moved to the side of the strip, then use the observed y-independence to extend our results away from the 
boundary. This is accomplished for both wired (fixed) and open (free) boundary conditions on the vertical ends. 

For the rectangle, we solve the density at the wall case with wired boundary conditions by use of vertex operator 
techniques. This derivation leads to expressions for the ratios of the various densities, including closed forms for the 
densities themselves on the boundary. 

In the limiting case of the semi-infinite strip we push the conformal calculation a bit further and obtain explicit 
formulas for all the various probability densities without the assumption of ^-independence, verifying the results just 
mentioned. This is accomplished for both wired (fixed) and open (free) boundary conditions on the vertical ends. 

All of our theoretical results agree very well with the numerical simulations. Neither finite-size effects nor correction- 
to-scaling operators make a significant numerical contribution. This was also observed in related work [6j], which 
includes comments on this situation. 

Taking the limit of a half-infinite strip, and then transforming the left-hand side to a semi-circle, gives a prediction 
for the behavior around a point as studied in [j| , where the radius e of the semi-circle is of the order of the lattice 
spacing. Thus, using our conformal field theory results, we are able to understand the finite-size effects for that 
problem, finding for example that the asymptotic behavior Co is approached as (e/r) 2 where r is the radial distance 
to the point. 

With open boundary conditions on the left- and right-hand sides of the rectangle, the behavior of C{z) is more 
complicated, and the factorization breaks down near the vertical sides. However, away from the sides, the function 
still approaches a constant Co, so there is factorization. Finally, with periodic boundary conditions on the horizontal 
sides and wired boundary conditions on the vertical ends, our simulations show approximate factorization everywhere, 
but a different asymptotic value for the constant: C\ » 1.022. 
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In this appendix, we verify that the solutions for the densities in section llll Bl indeed satisfy the null-state conditions 
imposed by the <\>\^. bcc operators. This proves that the y-independence is exact. 
Now the correlation functions representing the quantities of interest are of the form 
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APPENDIX: FULL DERIVATION OF SEMI-INFINITE STRIP DENSITIES. 



(01,2(O)0l, 2 (i)01/2 : o( z , Z)01,3(OO))S . 



(A.l) 



evaluated in the semi-infinite strip §. 

Consder the analogous quantity in the half-plane H: 



(01, 2 (1*1)01, 2("2)01/2, o{ w , w)01,3(«3))h . 



(A.2) 
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Conformal invariance implies that this correlation function may be written 

, -\~5/4% ( U2-U1 \ 1/3 ( (w - mi)(m 3 - u 2 ) (w - u^jus - u 2 ) \ 

[w-w) 1 F r, — , (A.3) 

\{u 3 - ui)(u 3 - u 2 ) J \{u 2 ~ ui)(u 3 - w) {u 2 —ui){u 3 —w)/ 

where F(rj, fj) is annihilated by the null state differential operators 

5/48 5/48 2/3 2d w 2d m 2d U2 2d U3 



3d Ul 2 and (A.4) 



(w — U\) 2 (w — Ml) 2 (1L3 — Ui) 2 w — Ul W — Ui U 2 — U\ U 3 — U\ 

Zd u 2 , (A.5) 



5/48 5/48 2/3 2d w 2d„ 2d Ul 2d U3 ., 



(W - U 2 ) 2 (w - U 2 ) 2 (U 3 - U 2 ) 2 W — Ui W - 1t 2 1*2 - Ui u 3 -u 2 

that arise from the </>i 2 operator at u\ and u 2 respectively. 

Next apply these two differential operators to (|A.3|1 and take the limit {u±, u 2 ,u 3 } — > {0, f,oo}. The correlation 
vanishes as u 3 goes to infinity. We account for this by replacing the second factor in (|A.3|) with the corresponding 
three-point function. Thus 

(&,2 (0)^,2(1)^/2,0 Kw)^i,3 (00)) n = fc(0)^i,2(l)feH)( w -i»)- 5 / 48 FKw) , (A.6) 

with F(w,w) satisfying 

Q = 5(w-w) 2 F _ 2{1 -wf dF _ 2(1 -w) 2 OF _ _ 2 d*F 
A8w 2 w 2 w dw w dw dw 2 

-w-w-v^-w-^m and ' 

5(w-w) 2 „ 2w 2 dF 2w 2 dF n „d 2 F „ d 2 F o „d 2 F 

= , gn V ' .. 2 F + - j- + - tt^ -3w 2 —- 6ww—— - 3w 2 — . (A.8) 

48(1 — iv) z (l — wy 1 — wow 1 — w dw dw z awow ow z 

Notice that these equations are interchanged by the transformation w — ► 1 — w, w — ► 1- w, i.e., there is mirror 
symmetry about the line u = 1/2. 

We transform these results from complex cooridinates in H into Cartesian coordinates in § using the mapping 

cosh(7T2:) + 1 

w(z) = , (A.9) 

so that 

(0 1 ,2(O)^,2(l)^/ 2 ,o(^^)0i,3(oo))s = (|!:§j) 5/96 M*) " w(z)T 5/4S F(w(z),w(z)) 

^M^) 2 + sH^y) 2 \ 5/96 H{Xiyh (A10) 



sinh(7ra;) 2 sin(7rw) 2 
where H(x, y) := F (w(x + iy), w(x — iy)). Using the relations 



d 2 d _ 1 r d .d 

dw tt sinh(7rz) dz tt sinh(7r(a; + iy)) \dx dy 

d 2d 1 / 8 . d , 

i— (A.12) 



dw 7rsinh(7rz) dz tt sinh(7r(a; — iy)) \dx dy 
we transform equation (|A.7j) for F(w,w) into one for H(x,y). Simplifying yields 

3sinh(7rx) 2 d 2 H 6 sinh(7nr) sin(7rw) d 2 H 



= 



7r 2 (cos(7ry) + cosh(7nr)) 2 9a; 2 7r 2 (cos(7ry) + cosh(7rx)) 2 dxdy 

3sin(7rw) 2 d 2 H (2 + cosh(7ra)(3 cos(7rw) + cosh(7rcc))) sinh(7ra) dH 

it 2 (cos(7ry) + cosh(7Tx)) 2 dy 2 7r(cos(7rw) + cosh(7re)) 3 dx 

(2 + cos(7rw)(cos(7r?/) + 3 cosh(7rx))) sin(7ry) dH 5 sin(7ry) 2 sinh(7rx) 2 



(A.13) 



7r(cos(7rw) + cosh(7ra;)) 3 dy 3(cos(7ry) + cosh(7ra;))' 



-H . 
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The other equation implied by (|A.8[) follows from this one by the mirror symmetry y — > 1 — y. 

Now the numerical evidence shows the three densities Pl(z),Pr(z) and Plr(z), that we calculate all take on the 
form 

Pi(z) = Qi(x)f(x,y), (A.14) 
with common function f(x,y). The function Pl(z) (see (l2~5j) and Fig|2~lb) may be written 

/„% , . . , .. /sinh(7ra;) 2 + sin(7rw) 2 \ 5 ^ 96 / sin(7ry) 2 \ ^ 6 , . „ „. 

(0i,2(O)0i, 2 (*)0i/ a ,oCM) s - ■ w \ 2 . / § ■ w J , . 7 ^ 2 (A.15) 

\ smn(7ra;) z sm(7ry) z y \smn(7nrj z + sm(7r?/J z / 

1 / sinh(7ra;) 2 sin(7ry) 2 N U 



simi^x) 1 / 3 \sinh(7ra;) 2 + sin(7ry) i 

This suggests setting 



11/96 

/(*. y) = ( iz„z ::;r;::., 2 ) ■ ( A - ie ) 



sinh(7ra;) 2 sin(7ry) z 
sinh(7ra;) 2 + sin(7n/) ; 



Comparing this to (|A.10[) further suggests that we let 

. . / sinh(7ra;) 2 sin(7rt/) 2 \ T , . 

= . , , J — ■ / \, (A.17) 

v yy \sinh(7rx) 2 +sm(7ry) 2 J v yy ^ ; 

It follows that if Pl(z), Pr{z) and Plr(z) do have a common y— dependence, then there must be solutions for J(x, y) 
that are independent of y and can be identified with Qr(x) and Qlr(x). We now demonstrate that this is indeed the 
case. 

Some algebra now shows that 

d 2 J d I d 2 1 d I 

= -18sinh 2 (7rx)— +37rsinh(27rx)— - 18sin 2 (7ry)— +37rsin(27ry)— + 10tt 2 J , (A.18) 
ox 2 ox oy z ay 

- as ' (A ' 19) 

These equations are surprisingly simple, and have definite parity under y — > 1 — y. 
Now (IA.19[) implies a solution 

J{x,y) = g x {x) + g 2 {y) , (A.20) 

an encouraging result since we expect solutions of precisely this form. 
We insert this into (|A.18[) via 

J(x, y) = g 3 {e~ 2 ™) + g 4 (sin(^)) (A.21) 

and express the result in terms of the new variables X :— e~ 27TX and Y :— sin(7ry). This leads to 

= (l0g 3 (X)-3(l-X)(7-5X)g 3 '(X)-l8X(l-X) 2 g 3 "(X)) (A.22) 
+ (10.94(F) + 6Y(1 + 2Y 2 ) 9i '{Y) - 18y 2 (l - Y 2 ) 9i "(Y)) . 

Because X and Y are independent, the two terms in (|A.22[) must either equal constants that sum to zero (which 
is the trivial solution since in that case J(x, y) = 0) or they must be independently equal to zero. This means that 
there are four linear solutions for J(x,y) that can be determined from two second order equations. 

The two solutions to the g 3 (X) equation are 

Qr(z) = U h smh(TTxy 1/3 e 27rx/3 2 F 1 (-2/3, -1/2, 5/6, e^ 2 ™) and (A.23) 
Qlr(z) = U h sin]i(7rx)- 1 / 3 e^/ 3 2 F 1 (-l/2,-l/3,7/6,e- 2 ^) , (A.24) 

which reproduce our expressions for Pr(z) and Plr(z) when combined with the prefactor f(x,y). For completeness 
we note that, as mentioned, Ql(z) — sinh(7rx) -1 / 3 with our particular choice of f(x,y). This completes the CFT 
proof of our observation that ratios of Pl(z), Pr(z) and Plr(z) are independent of y. 
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FIG. 21: Configurations in the semi-infinite strip: (a) Plr(z), (b) Pi R (z), (c) P L n(z), (d) Pf B (z), (e) P T g(z), (/) Ptb(z) 



Now this set of conformal operators also yields an additional set of correlation functions. Logically we expect that 
these solutions might be some of the corresponding quantities from the analogous limit of vertical crossing in an 
infinitcsimally short rectangle. The two solutions to the 34 (K) equation are given by 



Q TB (z) = n ft sin 5 / 3 ( 7 ry) 2 F 1 (4/3, 5/3, 2, sin 2 (^)) and 



Q iB {z) = n /l sm 5/3 (7rj/) 2 F 1 (4/3, 5/3,2, cos ; 



" 2 (f 



(A.25) 
(A.26) 



which represent configurations with the top and bottom conditioned not to belong to the same cluster, and the bulk 
point correlated with the top (Fig. [2Tb ) and bottom (Fig. |2"TH) sides respectively as in Fig. |2~T1 

These associations are based on leading term behavior as y — > and 1. For example, P T g(z) ~ y 2 and (1 — y) . 
The weight corresponds to the identity operator, which occurs in the bulk boundary OPE when the bulk operator 
approaches a fixed interval with which it is correlated. The weight 2 = h\ t 5 corresponds to the Fourtuin-Kastelyn 
four-leg operator, which we expect to appear when the bulk operator approaches a fixed boundary of different spin. 
This is best understood as the bulk operator pinching a dual cluster between itself and the boundary, a necessary 
condition in order that the boundary and bulk point have different spins. The inner and outer side of the dual 
cluster emanate from the boundary on both sides of the bulk point leading to a total of four Fourtuin-Kastelyn hulls. 
The combination of these two limits uniquely associates (|A.25I) with Fig. l2"T1 e while an identical argument fixes the 
association of (|A.26[) with Fig. |2TR 

For completeness we include 



Q TB {z) = sin^y)- 1 / 3 



(A.27) 



which follows from the correlation function 2(i)0i, 2(0)</>i/2, o( z : ^)) which was evaluated in 

We used (jA.23HA.24)) to construct the expression (|2~lj) for C(z) in the main text, but as we note in the introduction, 
the equivalent vertical quantity in this limit becomes 



C v (z) 



Ptb(z) 



^P t (z)P b (z)U v 



1 



(A.28) 



because II„ = 1 while Pt = Ptb + Ptb an d Pb = Ptb + Ptb both equal Ptb since II/j = 0. This is what we expect 
for crossing of a narrow rectangle. 
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